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ISOTHERMAL  DIFFUSION  IN  THE  DILUTE  RANGE  OF  THE 


SYSTEM  Ca  (OH)  2 — H3P04— H,,  O.  I - THEORY 
E.  C.  Moreno,  P.  R.  Patel  and  W.  E.  Brown 

Abstract 

The  equations  describing  isothermal  diffusion  in  the 
dilute  range  of  the  ternary  system,  containing  a weak  elec- 
trolyte, Ca (OH) 2-H3P04-H2 0,  are  derived  in  two  ways,  first, 
assuming  that  the  components  are  electroneutral  species  and, 
second,  considering  the  actual  ionic  species  present  in 
solution.  It  is  shown  that  the  two  models  are  thermodynamically 
equivalent.  The  theory  permits  the  calculation  of  the  four 
fundamental  diffusion  coefficients  (phenomenological  coeffi- 
cients) in  the  concentration  range  where  the  Debye-Htickel 
theory  suffices  for  the  calculation  of  ionic  activity  coeffi- 
cients. Therefore,  the  equations  can  be  used  to  test  the 
Onsager  reciprocal  relations  for  the  diffusion  process  in 
the  above  system.  The  ionic  model  allowed  derivation  of  ex- 
pressions for  calculating  practical  diffusion  coefficients  for 
the  electroneutral  components  from  the  limiting  equivalent 
conductances  of  the  ions  in  solution;  thus,  it  becomes  possible 
to  compare  experimental  diffusion  coefficients  with  those 
anticipated  theoretically. 
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relations,  it  is  necessary  to  calculate  the  partial 
derivatives  of  activity  coefficients  with  respect  to 
the  concentration  of  the  components.  Previous  treat- 
ments of  this  problem  in  ternary  systems  have  either 
assumed  invariance  of  activity  coefficients  with  con- 
centration10 or  have  used  assumptions  applicable  to 
symmetric  electrolytes11.  The  present  treatment  applies 
to  the  ternary  system  Ca (OH) 2 -H3P04-H2 0 in  the  concen- 
tration range  where  the  Debye-Huckel  theory  has  been 
shown1  2 9 1 3 to  describe  adequately  ionic  activity 
coefficients.  Derivations  will  be  given  for  two  models 
in  which  electroneutral  components  and  ionic  consti- 
tuents are  considered,  respectively.  It  will  be  shown 
that  these  two  models  are  thermodynamically  equivalent. 

2.0.  Procedure 
2.1.  Selection  of  model 

All  the  equilibrium  properties  of  the  ternary  system 
referred  to  in  this  paper  can  be  defined  in  terms  of  the 
components  Ca(0H)2/  H3P04/  and  H20.  For  diffusion 
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purposes,  the  number  of  components,  C1,  may  be  obtained 
from  the  expression  C'  = n-r,  in  which  n is  the  total 
number  of  chemical  species,  ionic  and  electroneutral, 
present  in  the  system  and  r is  the  number  of  restrictions 
operating  on  those  species.  Usually,  for  the  case  of 
strong  electrolytes  the  only  restriction  is  electro- 
neutrality. If  Ca(0H)2,  H3PO4,  and  H2 0 are  selected 
as  components,  in  general,  n would  comprise  eight  con- 
stituents: H20,  H3P04,  H2P07,  HPO",  P04,  H+ , Ca++,  and 

OH  . The  system  would  have  five  restrictions:  the 

three  expressions  for  the  three  ionization  constants  of 
H3PO4,  the  expression  for  the  ionization  constant  of 
H20,  and  the  condition  of  electroneutrality.  Therefore, 
the  system  is  ternary. 

It  is  neither  fruitful  nor  realistic  to  consider 
all  the  foregoing  constituents  and  restrictions 0 The 
derivations  would  become  unnecessarily  complicated? 
besides,  because  of  the  magnitudes  of  the  ionization 
constants  of  H3P04,  in  any  given  system  only  two  phos- 
phate containing  species  are  present  in  significant 


% 
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concentrations . Furthermore,  as  done  by  other  investi- 
gators10'1 4 , it  will  be  assumed  that  the  solvent,  H.,0, 
is  undissociated.  If  now,  neutral  components  that  do 
not  contain  the  ion  OH-  are  selected,  n will  be  equal 
to  five  (H20,  two  phosphate-containing  species,  Ca++ , 
and  H* ) , and  r will  be  two  (equilibrium  between  the  two 
phosphate-containing  species  and  electroneutrality) ; 
hence,  the  simplified  model  is  ternary  for  diffusion 
purposes.  The  two  phosphate  species  considered  here  are 
H3P04  and  PIP04  which  are  the  only  species  present  in 
significant  concentrations  in  the  experimental  pH  range 
of  4 to  6.5.  The  other  two  components,  besides  H20, 
correspond  to  electroneutral  combinations  of  the  ions 
considered  and  their  selection  is  arbitraryc  For  con- 
venience, CaHP04  and  H3P04  are  selected  here,  but  the 
treatment  would  apply  equally  well  to,  for  example, 
CaHP04  and  Ca(H2P04)2# 
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2.2  Derivations 


The  diffusion  coefficients  in  the  present  investi- 
gation were  obtained  in  a cell-fixed  frame  of  reference 
(a  porous  diaphragm  cell) . However,  for  the  range  of  the 
dilute  solutions  considered  here,  it  is  reasonable  to 
assume  that  the  partial  molar  volumes  of  the  components 
are  independent  of  concentration,.  Under  these  conditions, 
these  practical  diffusion  coefficients,  (Dt,)c,  become 
identical  with  those  referred  to  a constant-volume 
frame  of  reference,  (Dti)v  #l0#15  Furthermore,  the  solvent 
fixed  diffusion  coefficients,  (D13)0,  are  related  to 
the  volume-fixed  frame  of  reference10’14  by 


in  which  Cj  and  C0  are  the  concentrations  of  the  ith 
components  and  the  solvent,  respectively;  v0  and  vk , 
the  partial  molar  volumes  of  the  solvent  and  the  component 
K,  respectively;  the  summation  is  taken  over  q compo- 
nents, excluding  the  solvent.  It  is  apparent  that  at 


1 
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infinite  dilution,  (Dt , ) 0 = ; even  for  concentra- 

tions about  1M,  the  two  sets  of  coefficients  do  not 
differ  by  more  than  1 to  5 per  cent1  n . 

For  the  foregoing  reasons,  it  is  assumed  that  the 
derivations  shown  here,  although  strictly  applicable 
only  to  a solvent-fixed  frame  of  reference,  are  valid 
for  the  systems  under  consideration,  since  their  highest 
concentration  is  in  the  order  of  10~3M. 

Electroneutral  Model 

According  to  the  principles  of  non-equilibrium 
thermodynamics1 7 » 1 8 , the  flow  for  the  component  i,  Jl7 
in  a ternary  system  may  be  described  by 

j<  = - V Lj , i = o,  i,  2 (i) 

1 Lj 

3 =0 

in  which  , are  phenomenological  coefficients  and  the 
quantity  within  parenthesis  is  the  gradient  of  the 
chemical  potential  for  component  j along  x 7 the 
summation  in  (1)  includes  the  solvent  designated  by 


the  subscript  zero 
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From  the  Gibbs-Duhem  equation,  at  constant  tempera- 
ture and  pressure,  it  is  obtained 
a 

y CjdHj/TOC  = O (2) 

L.j 

3 ssO 


solution  of  equation  (2)  for  ay^/ax  and  substitution 

into  equations  (1)  gives 

2 

Jj  St  mm  \ X*!  j t = 1#  2 (3) 

3 = 1 

In  equations  (3),  the  flow  of  the  component  i is  given 
in  terms  of  the  chemical  potential  gradients  of  the 
solutes;  therefore,  the  phenomenological  coefficients, 
usually  called  "fundamental  diffusion  coefficients" 

Lt  3 , display  Onsager ' s reciprocal  relations1  4 , 

The  flow  Jt  in  equations  (3)  may  now  be  given  in 
terms  of  practical  diffusion  coefficients,  D1 s . In 
a ternary  system,  only  the  concentrations,  Cx  and  C2 , 
of  two  of  the  components  are  independent  variables,. 
Therefore, 


2 


3 = 0,  I,  2 


(4) 
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Substitution  of  equation  (4)  into  equation  (3)  gives 


Ji  = 


L>  j l ^ 


Sx 


(5) 


5 =1  k =1 


or 


Ji 


2 

r1 

“ ) ^ i k 

k =1 


(6) 


in  which  the  practical  diffusion  coefficients,  Dik  are 
given  by 


D 


i k 


3 =1 


BCk 


1 =1»  2 
k =1,  2 


,(7) 


The  four  practical  diffusion  coefficients  Dlkfcan  be 
obtained  through  a minimum  of  two  diffusion  experiments, 
by  the  use  of  equations  (6)  with  the  appropriate  boundary 
conditions.  From  equations  (7) , the  following  relations 
may  be  obtained  by  solving  for  the  Lt , : 


^11  (^1  1 P'2  2 ^12^1  ^/S 

^12  = (^1  2 M*1  1 " D1  1 1^1  2 ) /S 

(8) 

^21  ~ ^2  1^2  2 “ ^2  2 M»2  1 ) /S 
■^2  2 — (^2  2^11  " ^2  1 ^12^/S 

in  which  ^ i stand  for  the  partial  derivatives 
and  the  denominator  s is  given  by 


s 


M*l  l po  p. 


2 1 


Equations  (8),  previously  used  by  other  investigators1  4 * 1 
permit  to  calculate  the  fundamental  diffusion  coefficients, 
L. j , from  practical  diffusion  coefficients  that  are 
determined  experimentally,  provided  that  the  values  of 
the  various  partial  derivatives  are  known.  The  method 
for  evaluation  of  these  derivatives  is  given  next. 

The  chemical  potential  of  the  ith  component,  p,t , is 
defined  in  terms  of  its  activity  al , by 
M*i  - M*?  + RT  In  at 

in  which  p,  J is  a constant  under  constant  temperature,  T, 
and  pressure,  and  R is  the  gas  constant.  The  partial 
derivatives,  i , are  given  by 


RT 


i = 1,  2 

3 = 1*  2 


(9) 


in  which  at  , represents  the  partial  derivatives  c^/dC,  . 
Thus,  the  evaluation  of  p,t  , is  accomplished  if  appropriate 
expressions  of  a4 3 are  available. 
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As  shown  in  equations  (46)  and  (47)  in  Appendix 
lf  the  activities,  ax  and  a2  , for  components  1 and  2 
(CaHP04  and  H3P04,  respectively)  are  functions  of  their 
concentrations  Cx  and  C2 . Other  parameters  in  those 
expressions  (hydrogen  activity  and  ionic  activity 
coefficients)  are  also  functions  of  Cx  and  C2  0 There- 
fore 


3 £ k 


in  which  H represents  the  ionic  activity  of  hydrogen,  and 
Y^  represents  the  activity  coefficient  of  each  of  the  i 
ions  assumed  to  be  present  in  the  system,  Ca++,  H*  , 

H2P04/  and  HP04. 

The  electroneutrality  function  E,  equation  (42) 
in  Appendix  1,  which  will  be  used  later,  is  defined  on 
the  basis  that  only  these  ions  are  present.  (The 
presence  of  ion  pairs  is  ignored.)  In  the  subsequent 
equations  through  equation  (19) , indexes  1 / j , and  k , 
assume  the  same  values  as  in  equation  (10) . 


11- 


In  the  dilute  range  of  the  ternary  system  under 
consideration,  it  has  been  shown12’13  that  the  Debye- 
Hiickel  theory  provides  an  adequate  expression  for  the 
calculation  of  ionic  activity  coefficients.  Thus, 
y ^ are  explicit  functions  of  the  ionic  strength  I. 
Therefore,  the  second  factor  in  the  summation  term  of 
equation  (10)  becomes 


By  substitution  of  equation  (11)  into  equation  (10)  and 
making 


(12) 
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By  differentiation  of  equation  (43)  in  Appendix  1,  it 


follows  that 


(13) 


By  substitution  of  equation  (11)  into  (13)  and  making 


equation  (13)  becomes 


vSC  ,yc. 

v u k 


(14) 


Solving  equation  (14)  for  (dl/dC,)c, 

"dl 


BC 


IVdl 

['vsc, 


y|V  H- 


+ r— " i 

V3H^'C  ,,  v Vac.y. J 


(l-I*)"1  (15) 
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Substitution  of  equation  (15)  into  equation  (12)  gives 
an  expression  for  at 3 in  terms  of  partial  derivatives 
that  can  be  readily  calculated  from  experimental 


measurements,  except  for  the 
derivation  for  (?*H/?>C,)r  is 

J k 

Differentiation  of  the  < 
E,  equation  (42)  in  Appendix 


By  substitution  of  equation 
and  making 


derivative  ( ) c . The 

J k 

given  next. 

lectroneutrality  function 
1,  gives 


(16) 


(11)  into  equation  (16) 
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equation  (16)  becomes 


Substitution  of  equation  (15)  into  equation  (17)  gives 
the  solution  for  (dH/6Cj)Ck, 


Substitution  of  equation  (18)  into  equation  (15)  gives 


The  right-hand  terms  in  equations  (18)  and  (19)  involve 
quantities  that  can  be  readily  calculated  from  experi- 
mental measurements.  The  actual  expressions  used  for 
such  calculations  are  given  in  Appendix  2.  Substitution  of 
equations  (11) , (18)  and  (19)  into  equation  (10)  gives 

the  needed  expression  for  the  evaluation  of  at , ; then, 
by  the  use  of  equation  (9),  evaluation  of  u is 
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accomplished , and  substitution  of  the  (jq  into  equations 
(8)  permits  the  calculation  of  the  fundamental  diffusion 
coefficients,  , and,  therefore,  the  testing  of  the 

Onsager ' s reciprocal  relations, 

Ionic  Model 

In  the  foregoing  treatment,  the  components  have 
been  considered  as  neutral  species,  although  in  defin- 
ing their  activities,  their  ionic  nature  was  taken  into 
account.  Consideration  of  an  ionic  model  for  the  present 
system  is  advantageous  in  that  practical  diffusion 
coefficients  can  be  calculated  for  the  electroneutral 
components  on  the  basis  of  ionic  conductances1  0 . It 
is  then  possible  to  compare  the  practical  diffusion 
coefficients  derived  here  with  those  obtained  experi- 


mentally 
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The  ionic  model  selected  here  is  given  in  equations 
(20)  to  (22) 


CaHP04  - Ca++  4-  HP04 

(20) 

H3P04  ->  2H+  4-  HP04 

(21) 

H2P04  ^ if  4-  HPOj 

(22) 

This  ionic  model  implies  that  there  is  no  undissociated 
H3P04/  an  assumption  which  is  valid  for  the  pH  range 
considered  hereQ  It  will  be  shown  now  that  this  system 
is  thermodynamically  equivalent  to  the  one  in  which 
CaHP04  and  H3P04  are  considered  as  the  neutral  components. 
In  the  following  derivations,  i represents  the  consti- 
tuents CaHP04 , H3P04/  and  H2P04  and  j represents  the 
ionic  species  in  equations  (20)  to  (22) . The  numerals 
1,  2 and  3 are  used  for  these  three  constituents,  respec- 
tively and  the  numerals  4,  5 and  6 represent  the  ions 
Ca++ , HP04  and  H+ , respectively. 

The  stoichiometric  coefficients  in  the  right-hand 
terms  of  equations  (20)  to  (22)  are  represented  by  V{ , . 

For  example,  4 is  the  stoichiometric  coefficient  for 
the  Ca++  that  comes  from  CaHP04  and  vx 5 is  for  the 
HP04  that  comes  from  the  CaHP04#  The  numerical  values 
for  the  's  are:  v1 4 = vis  = v2 B = v3s  = v36  = 1; 

v16  = 


V24 


v3  4 — 0 ; v2  6 


2. 


-17- 


On  the  basis  of  the  ionic  model,  the  entropy 
production,  <7 T,  in  the  diffusion  system  is  given  by1  ' ’ 2 


in  which  Ji  represents  the  diffusion  flow  of  the  jth  ion 
and  its  generalized  driving  force,  in  this  case,  its 
electrochemical  potential  gradient;  cr  is  the  rate  of 
entropy  production  at  temperature  T. 

We  consider  now  as  new  generalized  forces  the 
chemical  potential  gradients,  X4 , of  the  neutral  com- 
ponents 1 and  2 and  the  electrochemical  potential  X3 0 
The  transformation  equations  for  the  new  forces  are 


(23) 


6 


Xj  — \ Vij  X 


u 


( I =1  , 2) 


(24) 


6 


3 “4 
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which  are  linear  combinations  of  the  old  forces  Xj . 

The  transformation  equations  for  the  new  flows,  derived 
from  the  mass  balance  condition  in  equations  (20)  to 
(22) , are 


J 


3 


^3 


l J=4  tO  6 ) 


1 =1 


(2  5) 


Substitution  of  equation  (25)  into  (23)  gives 


But,  according  to  equations  (24)  the  summation  term 
within  parenthesis  equals  X3 , therefore,  the  quantity 
within  parenthesis  vanishes  and 

2 

<jT  = V J,  Xj  (27) 

i =1 

hence,  the  ionic  model  is  thermodynamically  equivalent 
to  the  diffusion  system  in  which  the  neutral  compounds 
CaHP04  and  H3P04  are  considered  as  components. 
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The  procedure  to  calculate  fundamental  diffusion 
coefficients  for  neutral  components  from  ionic  conduct- 
ances has  been  given  in  detail  by  Wendt1  0 and  it  was 
tested  by  the  same  author  with  diffusion  data21  for 
the  system  H2 S04-Na2 S04-H20o  For  this  reason,  only  a 
summary  of  our  derivation  is  given  next  in  which  we 
apply  Wendt's  treatment  to  the  present  system.  We  use 
the  $me  terminology  except  that  his  Lt,  (for  ions) 
corresponds  to  our  Ft  3 and  his  a* i correspond  to  our  at . • 
The  diffusion  flows  J3  of  the  j ions  are  related 
to  the  chemical  potential  gradients,  Xn , of  the  two 
neutral  components  by 
2 

= V_1  l 

n.  =1  k =3 

in  which  v is  the  value  of  the  determinant  of  the  vt 2 
matrix  (equal  to  -1  in  the  present  case);  ak n , the 
elements  of  the  matrix  [ockn]/  hava  the  numerical  values 


1 

“ O3  2 “ Ct5i  — 

06  1 = 0/ 

O3  3 “ 041  “ Os  2 — 1, 

05  2 ~ 

a6  3 

■*-  / a4  3 2 , 

0,5  3 ="*2; 

the  coefficients  Qjk 

are 

given  by 


c 

akn  ^Jk  ( J =3  tO  6 ) (28) 
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- i j 3 I k k ak3 


O 


q — 3 


! =n  to  o 


(29) 

k =3  to  6 


In  equation  (29).,„  6,k  is  the  Kronecker  delta  (it  has 

the  value  1 for  j=k,  and  zero  for  jtk) , and  the  P's 

are  the  main  terms  of  the  fundamental  diffusion  coefficients 

for  the  subscripted  ions,  3 to  60  Equation  (29)  is 

consistent  with  Wendt's  assumption10  that  the  coefficients 

F3k  are  zero  for  j=j=k0  These  ionic  coefficients  can 

be  calculated  from 

rn  = 1°  c,|  z,|  f2  x io7  i=3  to  s (30) 

in  which  X°  is  the  limiting  equivalent  conductance  of 
the  ion  j,  C,  its  molarity,  ] Z,|  the  absolute  value  of 
its  valance  and  F is  the  Faraday  constant  (96,493 
coulombs  per  equivalent  of  charge).  Values  used  for  \ ° 

r -i  • 1 • V 

in  this  investigation  are  given  in  Appendix  1. 
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Expanding  equation  (25)  and  solving  for  Jx  and 
J2  (using  the  numerical  values  for  the  corresponding 


) we  get 

i = 

(31) 

2 ~ 2 ( J3  + dr  ) 

(32) 

The  expressions  for  the  ionic  flows  J3 , J4  and  J6 
given  by  equation  (28)  can  now  be  substituted  into  equa' 
tions  (31)  and  (32)  to  obtain 


2 6 

J"i  — “ ^ ^ ak  r_  Xn  Cl 4.  k 

n =1  k =3 


(32) 


2 6 


2 6 


J2  = 


^ ak  n ^6  k " 2 ^ ^ akn  ^3k 


(33) 


n =1  k =3 


n —1  k =3 


But  from  equation  (3)  we  can  also  write 


(34) 


(35) 


n =1 
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Through  a comparison  of  the  coefficients  of  Xn  in 
equations  (34)  and  (35)  with  those  in  equations  (32) 
and  (33)  respectively,  it  is  obtained 

6 

Lin  ^ &k  n ^'4k 

k =3 

n=l,  3 (36) 

6 

^3  n = i ^ ak  n ( ^6  k + ^3k  ) 

k =3 

Equations  (36)  together  with  (29)  and  (30),  permit  the 
calculation  of  theoretical  fundamental  diffusion 
coefficients,  Lt } , from  a knowledge  of  the  limiting  con- 
ductance of  the  ions  involved,  subject  to  the  assumption 
adopted  (ideal  ionic  behavior) . These  coefficients  are 
compared  with  those  calculated  from  experimental  data  in 
the  results  section;  likewise,  a comparison  is  made 
between  the  values  for  the  practical  diffusion  coefficients, 
Dj  3 , obtained  experimentally  and  those  calculated  by 
substituting  the  , obtained  in  equation  (36)  into 
equation  (8) , using  the  expressions  for  Hq , previously 


derived 


Symbols : 


APPENDIX  1 


C-l  = molar  concentration  of  CaHP04 

C2  = molar  concentration  of  H3P04  used  in  preparation 

of  solution 

P0  = molar  concentration  of  undissociated  H3P04 

P1  = molar  concentration  of  H2P04 

P2  = molar  concentration  of  HP04 

y1  = molar  activity  coefficient  for  H2P04 

y2  = molar  activity  coefficient  for  HPO^ 

y 3 = molar  activity  coefficient  for  Ca+'r 

y 4 = molar  activity  coefficient  for  H+ 

[H]  ~ molar  concentration  of  H+ 

H = activity  of  H+ 

Kx  = first  ionization  constant  of  H3P04/  7.10e  x 10~3  (22) 

K2  = second  ionization  constant  of  H3P04,  6.33e  x 
10-8  (23) 

(P2)  = activity  of  HP04 

\ ° 3 = limiting  equivalent  conductance  of  H2P04/  32.3  (24) 

\°4  = limiting  equivalent  conductance  of  Ca++ , 59.50  (25) 

\°5  = limiting  equivalent  conductance  of  HP04/  43.7  (26) 

X °s  = limiting  equivalent  conductance  of  H+ , 349. 8X  (25) 


APPENDIX  1 (continued  2) 


The  total  molar  concentration  of  phosphorus  in  the 
solutions  of  CaHP04  in  dilute  H3P04  is  given  by  (C]  f C2  ) . 
The  balance  equation  for  phosphorus,  for  solutions  with 
pH  values  below  7 (P04  is  negligible) , may  be  written 

(Cx  4-  C2 ) = P0  f P,  4-  P2  (37) 

in  which  P0  is  the  molar  concentration  of  undissociated 
H3PO4.  The  quantities  P0  and  PT  may  be  written  in  terms 
of  P2  by  the  use  of  the  expressions  for  the  first  and 
second  ionizations  of  H3P04 7 then,  solving  equation  (37) 
for  P2  gives 


“ (Cj.  4- 


C2) 


( 


4- 


H 


Kp  Yi 


4- 


H' 


Ki  K2 


(38) 


In  the  pH  range  4 to  5,  the  first  and  third  terms 
within  the  parenthesis  in  the  denominator  of  equation 
(38)  are  negligible  compared  to  the  second  one.  Between 
pH  5 and  7 only  the  first  and  second  terms  in  the  paren- 
thesis need  to  be  consideredo  The  maximum  error  intro- 
duced in  the  value  of  P2  by  these  simplifications  amounts 


to  1 per  cent 


APPENDIX  1 (continued  2) 


Calling  fs  the  quantity  within  parenthesis  in  the 
denominator  of  equation  (38) 


PE  = (Cj  -I-  C2)/vaf. 


(39) 


and  by  substitution  of  equation  (39)  into  the  expression 
for  the  second  ionization  of  H3P04 , it  is  obtained 


Pl  = (Cl  + C2)H/f27lK. 


(40) 


The  function  E,  that  defines  electroneutrality,  is 


E = 0 = 2Ci  + — - 2 P,  - P, 

Y4 


(41) 


substitution  of  equations  (39)  and  (40)  into  (41)  gives 


E = O = 2Ci  + 


H 


YjK. 


+ 


(42) 


The  ionic  strength,  I,  is  given  by 


H 


I - 2CX  + 2y  + 


Cx  + C j 
2 f o 


H 


+ 


ViKs  Vs 


(43) 


For  brevity  in  the  expressions  given  in  Appendix  2, 
we  shall  define  the  quantities  within  square  brackets 
in  equations  (42)  and  (43)  as 


APPENDIX  1 (continued  3) 


Q = 


H 

^2  Vl 


"f* 


V: 


(44) 


G = 


H 

K2Yi 


+ 


Ys 


(45) 


The  activities  of  CaHPCh  and  H3P04  in  solution 


are  given  by  equations  (46)  and  (47) , respectively. 


(Ci2  + C2 C-l  ) y3 


a,  = 


(46) 


a,  = 


(Ci  + c2)h 

fiK, 


(47) 


The  function  fj  in  equation  (47)  is  defined  by 


f =—  + 

Kj  HYe 


+ 


(48) 


Yi 


The  ionic  activity  coefficients, 
from  the  Debye-Huckel  theory 


i ' 


are  calculated 


AZ|  /I 

In  y.  = - t— 77 

H l+Ba  ,/i 


(49) 


in  which  the  constants  A and  B,  dependent  on  temperature 
and  the  dielectric  constant  of  water,  have  numerical 
values  of  0o5092  and  0.3286  x 10s , respectively  (27)  at 

25°C;  Z is  the  valence  of  the  j£th  ion,  and  a*  its 

H/ 


APPENDIX  1 (continued  4) 

distance  of  closest  approach.  Numerical  values  for 
(28)  are:  Ca+ h , 6 x 10“* cm;  H+ , 9 x 10“* cm;  HP04 

and  H2  P04  , 4 x 10“ft 


cm. 


APPENDIX  2 


Expressions  for  the  calculation  of  the 
derivatives  used  in  text. 

The  expressions  given  in  this  Appendix 
on  the  equations  in  Appendix  1. 


dY^  $ (I) 

dl  - yl  21 (1+Ba^/l 


in  which 


$ (I) 


AZ^  /I 
1+Ba^/X 


2 


+ 


G 

2fs 


1 (Ci  + C3) 

2yT  + 2f2Ka 


partial 
are  based 


(50) 


(51) 


(52) 


(53) 


APPENDIX  2 (continued  2) 


(Ci  4-  c3) 

2 f 2K2  y i 


■ 

4 


(54) 


(55) 


= ( 2CX  + Cs)  y3 

y i f2 


(56) 


Y3 

f« 


(57) 


■3  a, 

iiT 


3» 


(cx  + cx  Cs  ) y3  ,2H  + ]_N 
f|K2  % Yly 


(58) 


3 a, 


(C?  + CXCS) y3H 


C ^ H » 7 4= , 


K 


2 

2 Yl 


(59) 


(C|  + C* C2)y3 

y^2  flvl 


(60) 


APPENDIX  2 (continued  3) 


3*1 


d y 


(Cf  + CjCg) 


3 ■'  C 


I f c 2 f H » 


(61) 


H 

fuq 


(62) 


(Ci  + C2 ) ^ K2  H x 
ffKj  V1  + 5^  " K,J 


(63) 


^a2\ 


Y 

2’  ifl 


(Ci  + C2)H 

fi2K4Yi 


(64) 


(Ci  + C2)K2 

f?KiY| 


(65) 


(66) 


1» 


H,  Yjl 


(67) 


APPENDIX  2 (continued  4) 


Y* 


1 


Y 


4 


(Ci  + ca)  ^ _ Q_  - 

f2KaYi  v ” f,  y 


(68) 


(C,  + C2)H 


f2K2v; 


(69) 


(c^  + cs ) Q ^ 

f2Y,  V " f3- 


(70) 


r— ) 

V3y  J 


**  c 2»  H>  \<= 4 


H 

Y: 


(71) 
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